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We discuss the beiiaviour of the Dirac fermions in a general spherically symmetric black hole 
background with a non-trivial topology of the event horizon. Both massive and massless cases are 
taken into account. The analytical studies of intermediate and late-time behaviour of massive Dirac 
hair in the background of a black hole with a global monopole and dilaton black hole pierced by a 
cosmic string will be conducted. It was revealed that in the case of a global monopole swallowed by 
a static black hole the intermediate late-time behaviour depends on the mass of the Dirac field, the 
multiple number of the wave mode and the global monopole parameter. The late-time behaviour 
^ . is quite independent of these factors and has the decay rate proportional to t~^^^ . As far as the 

'^1^ ' black hole pierced by a cosmic string is concerned the intermediate late-time behaviour depends 

\ only on the hair mass and the multipole number of the wave mode while the late-time behaviour 

. dependence is the same as in the previous case. The main modification stems from the topology of 

the 5*^ sphere pierced by a cosmic string. This factor modifies the eigenvalues of the Dirac operator 
acting on the transverse manifold. 

PACS numbers: 04.50.-|-h 

^ ■ I. INTRODUCTION 

>. . 

CS) ^ During past decades the study of fermions behaviour in various backgrounds has attracted much interests. Exact 
solution of Dirac equation in curved spacetime is a very useful tool in investigations of physical properties of particles 
in adequate spacetimes. As far as black holes is concerned, they are in this category, because of the fact that their 
. better understanding acquires a detailed study of various matter fields in its vicinity (see [l| and references therein) . 
' Dirac fermion fields were studied in the context of Einstein- Yang-Mills background 2] found by Bartnik and McKinnon 
, 3 . Fermion fields were analyzed in the near horizon limit of an extreme Kerr black hole |^] . It turned out that the 
' extreme Reissner-Nordstrom (RN) case Q having both magnetic and electric charges was destroyed in the presence 
. .\ of a massless fermion field. It was also shown generalizing previous results fT*!, that the only black hole solution of 
^ ■ four-spinor Einstein-dilaton- Yang-Mills equations were those for which the spinors vanished identically outside black 
hole. It physically means that the Dirac particles either enter the black hole or escape to infinity. In Refs.fl, |^ Dirac 
^ fields were considered in Bertotti-Robinson spacetime. On the other hand, Dirac fields as a cosmological solution 
' with a homogeneous Yang-Mills fields acting as an energy source were analyzed in Ref . [lOj . Dirac equation, via the 
Newman-Penrose formalism, in the context of Robertson- Walker spacetime was considered in [ll|. Recently, the 
behaviour of massive Dirac fields on the background of a charged de Sitter black holes was elaborated in Ref. [l3 • 

Regardless of the details of the gravitational collapse, the resultant black hole is described only by few parameters 
such as mass, charge and angular momentum. It is the famous Wheeler dictum that black holes have no hair. But 
the tantalizing question is how the loss of black hole hair takes place. The first contributions towards this problem 
were given in Refs.fl^, where it was shown that the late-time behavior is dominated by the factor t-(2'+3)^ foj. 
each multipole moment I and the decay-rate along null infinity and along the future event horizon was governed by 
the power laws and where u and v were the outgoing Eddington-Finkelstein (ED) and ingoing ED 

coordinates. In Refs.[ll]-[l3] charged scalar hair decay cases were considered. It turns out that a charged scalar hair 
decayed slower than a neutral one. The nearly extreme RN spacetime was treated in Ref. where it was shown 
that the inverse power law behaviour of the dominant asymptotic tail is of the form i^^/^ sin(mt), being independent 
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of I. The asymptotic tail behaviour of a massive scalar field was also studied in Schwarzschild spacetime The 
oscillatory tail of the scalar field has the decay rate of t~^^^ at asymptotically late time. The power-law tails in the 
evolution of a c harg ed massless and self-interacting scalar field around a fixed background of dilaton black hole was 
studied in Refs.(2Ct [2]|. where the inverse power-law relaxation of the fields at future timelike infinity, future null 
infinity and along the outer horizon of the considered black hole was found. It was also envisaged numerically that at 
very late times the oscillatory tail decay is of the form t^^^^. These results were confirmed analytically in the theory 
with arbitrary coupling constant between dilaton and J7(l)-gauge fields [2^ . 

Much attention was also paid to the late-time behaviour of massive Dirac fields. They were studied in the spacetime 
of Schwarzschild black hole [2^ as well as in RN black hole background [l^l • The stationary axisymmetric black hole 
case was studied numerically in Ref. j25| and it was found that in the case of Kerr black hole the oscillatory inverse- 
power law of the dominant asymptotic tail behaviour was approximately depicted by the relation t~^/^ sm(mt). In 
Ref. [2^ the decay of a charged massive Dirac hair in the background of Kerr-Newmann black hole was investigated. 
It turned out that the intermediate late-time behaviour of the fields under consideration was dominated by an inverse 
power-law decaying tail without any oscillation. 

In Ref. ^27*1 the analytical studies concerning the intermediate and late-time decay pattern of massive Dirac hair on 
the dilaton black hole were conducted. Dilaton black hole constitutes a static spherically symmetric solution of the 
theory being the low-energy limit of the string theory with arbitrary coupling constant a. 

Massive vector field obeying the Proca equation of motion in the background of Schwarzschild black hole was studied 
in |2SJ. It was revealed that at intermediate late times, three functions characterizing the field have different decay 
law depending on the multipole number I. On the contrary, the late-time behaviour is independent on I, i.e., the 
late-time decay law is proportional to sm{mt). 

Recently, there has been renewed interest in brane models in which our Universe is represented as (3-fl)-dimensional 
submanifold living in higher dimensional spacetime as well as brane black holes. The decay of massive Dirac hair on 
a brane black hole was considered in [29|]. 

It might hap pened that at the beginning of our universe several phase transitions lead to the topological defect 
formations [30| • The interaction of topological defects such as global monopoles or cosmic strings with black holes has 
very interesting features. For instance, the black hole global monopole system has an unusual topological property of 
possessing a solid deficit angle. The physical characteristics of the above system were widely studied in literature |3l| . 
The decay of massive scalar hair in the background of a Schwarzschild black hole with global monopole was studied in 
Ref. [3^ . It happened that the topological defects makes the massive scalar field hair decay faster in the intermediate 
regime comparing to the decay of such hair on the black holes without defects. On the other hand, the late-time 
behaviour was unaffected by the presence of global monopole. The Schwarzschild black hole global monopole system 
coupled to scalar fields was elaborated in Ref. [33] . 

On the other hand, cosmic strings and cosmic string black hole systems also acquire much interest. Assuming 
a distributional mass source the metric of this system was derived in [s^ (the so-called thin string limit). Later, 
it was revealed that it constituted the limit of much more realistic situation when a black hole was pierced by a 
Nielsen-Olesen vortex [ssj. It turned out that for some range of black hole parameters extremal black hole expelled 
the vortex (the so called Meissner effect). Moreover, extremal black holes in dilaton gravity always expel Higgs fields 
from their interiors (36| . 

It will be not amiss to study a much more realistic case than the scalar fields, i.e., the behaviour of fermion fields on 
the background of black holes with non-trivial topologies of their event horizons caused, e.g., by topological defects 
and the influence of the topological defects' parameters on the intermediate and late-time decay of massive fermion 
hair. 

The main aim of our work is to elaborate the solution of the Dirac equation in the background of a general spherically 
symmetric black hole spacetime. We shall try to generalize in some way the attitude presented in 2] studying zero 
modes and A: > modes of massless and massive Dirac equation in topology non-trivial spacetimes. First, we shall 
take into account the case of a Dirac equation supplemented by the Yang-Mills term. We assume the complete 
separation of the degrees of freedom of the considered fields. In our work, we shall consider the case of a black hole 
global monopole system which implies 

, 2 / 2GM\ , 2 dr^ ,2 2 ,r.2 

ds^ = -[l \dt^ + -r ^ + b^r^dn^, (1) 

V ^ / / 1 2GM \ 



as well as the black hole pierced by a cosmic string, which line element in the thin string approximation of the 
Nielsen-Olesen vortex passing through the spherically symmetric static dilaton black hole. We shall consider the 
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dilaton gravity with an arbitrary coupling constant a. Namely, the line element in question yields 

ds' = -(l-':±)(l-—] ^dt" + — — + R\r)dn\ (2) 



r \ r 



1-^1(1-^ 



l + Q 



where B? [r) — { 1 — — \ and the transverse two-dimensional manifold is given by dVl^ ~ dO^ + s\v? 



The parameter B is related to the linear mass density of a cosmic string passing through the dilaton black hole. Of 
course, for dilaton black hole without a cosmic string one should have B = 1. r+ and r_ are the outer and the inner 
event horizons of the black hole. They are related to mass M and electric charge Q of the black hole in the following 
way: 

e-^==(l-^)^, 2M^., + i^._, Q^-^. (3) 

Then, we proceed to analytical discussion of the decay of massive fermion hair on the aforementioned backgrounds. 

The layout of our paper is as follows. In Sec. II we first briefly review the behaviour of Dirac fermions in the 
background of a general spherically symmetric black hole. We analyze the massless case, zero modes of Dirac equation 
in the near-horizon limit for both extremal and nonextremal black holes with non-trivial topology of the event horizon. 
Then, we take into account the Dirac fermion fields for fc > and show that the equations in question can be 
decoupled to the system of second order differential equations. In the next subsection we treat the massive case of 
Dirac fermions. It also happens that their equations decouple to the system of second order differential ones with the 
so-called supersymmetric potentials. In Sec. Ill we gave the analytic arguments concerning the decay of Dirac massive 
hair in the backgrounds of topology non-trivial black holes. Finally, we conclude our investigations in Sec. IV. 

II. FERMIONS IN THE GENERAL SPHERICAL BACKGROUND 

In this section we shall begin our analysis by considering the behaviour of Dirac fermion fields in a general spherically 
symmetric black hole background. The metric corresponding to the aforementioned spacetime will be given by the 
expression 

ds^ = -A{r)^dt^ + Birfdr"^ + C{rY dfl^ , (4) 

where all the metric functions have the r-coordinate dependence. The transverse metric d^l^ depends neither on 
r-coordinate nor on t. In this background, we refine our studies to the solution of the massless Dirac equation with 
Yang-Mills potential provided by 

ir{\/^~tXH^)^ = 0, (5) 

where is the covariant derivative, =9^-1- \/2ijj^°'^"jalb and oj^""^ are the connection one-forms. The Dirac 
matrices satisfy {7", 7^} = —2ri°'^. They are given in the explicit form as follows: 

where / stands for the identity matrix while cr'^ are the Pauli spin matrices. The vector potential of the Yang-Mills 
type appearing in Eq.([S]) will be chosen as 

H. = ^e,,, r^ (7) 
where A is the coupling constant, r'^ is the generator of the group SU{2) and the unit normal vector. The spinor 
-0 can be decomposed into its left and right chiral component, i.e., ip — ( ^''^ ) . It is well known (see, e.g., [sJl) that 



for massless case the two chiralities decouple. It enables one to rewrite the underlying Dirac equations in the form as 



dt^R - tX A{r) a^'H^ ^jr + ^ P^jr + ^ -A-i_a,(C(r) A(r) 2 = 0, (8) 



Ml + ^X A{r) a^H, - ^ - ^ ^!ila,.(C(r) A(r)^ 4,^) = 0, (9) 
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where p is the Dirac operator on the transverse manifold. 



A. Zero modes of the Dirac equation 

First of all we shall consider the s-wave case. Because of the fact that we are looking for the massless spinor 
solution we can restrict our attention to the case when ij) = -0^, without loss of generality. In our considerations we 
shall use the hedgehog spinor ansatz (see e.g., (ssl [39j) in which the spinor function ip will be spanned by two states 

^ iTf (^o) (l) (o) ^'^^ ^ '^^^'^ ^o^'^O'^^'"' '^il^ ^^^^ property (crfe + Tk)xi = 0. 

Consequently, the considered Dirac spinor ■0 may be written in the form as 

^ = C(r)-iA(r)-V2 f(^t,r) xi + C{r)-^ A{r)-^'^ 9{t,r) X2. (10) 

Properties of the hedgehog spinors allow one ^ to find that the transverse Dirac operator will act as^ xi — ~X2 and 
P X2 — Xi- Moreover one finds that na x f xi = ^'^iX2 and na x f X2 = 2jxi- The above properties of the hedgehog 
spinors help us to rewrite the Dirac equation for the s-wave sector in the form 

(-■;";) ^c-)) ,,o. (12, 

r C{r) 

where for a convenience we have introduced the tortoise coordinate connected with the r-coordinate by the relation 
dr^,/dr = B{r) / A{r). By virtue of Eqs.(fTT |) - ([T2l) we can readily conclude that they admit a zero energy bound state 
given by 

It will be interesting to find the near-horizon behaviour of the Dirac fermion fields. In the case of a nonextremal black 
hole one can expand the metric coefficients A{r) and B{r) in the vicinity of the black hole event horizon. They will 
be provided by the following expressions; 

A{rf ~ A'{r+){r~r+), B{rf ~ B' {r+){r - r+y^ . (14) 

Then, making a change of variables given by the relations 



P'-^B'ir,)ir-r,l T=\^^^t^ (15) 

it can be shown that the line element describing the near horizon geometry of the nonextremal black hole can be cast 
in the form as 

= ~p2^r2 + dp' + C(r+)2(d02 ^ gij^2 Q^^2^_ (^g) 

Thus, one can approximate the spacetime in the vicinity of the nonextremal black hole event horizon by the Rindler 
line element. On this account, having in mind Eas.(|10p and (|13p. one can verify that the spinor function can be 
approximated by the following expression: 

V'- ^-1 r- (17) 

C{r+) A'{r+)i (r - r+)i 

One concludes that it behaves as (r — r+)^^/* near the event horizon r+. If we shall consider the global monopole 
black hole spacetime then C(r+) — br^^ where h is the global monopole parameter. Just, the bigger b we have the 
more divergent is the spinor wave function. In the case of a dilaton black hole, one can remark that the a-coupling 
constant will also trigger the divergence of the spinor function near the black hole event horizon. 

On the other hand, in the extreme black hole case, when the outer black hole event horizon is equal to the inner 
one, r+ = r_, the metric functions imply the following: 

A{rf ^ \A"{r+){r - r+)\ B{rf ^ ii?"(r+)(r - r+r\ (18) 
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Let us use a suitable change of the coordinates given by 



_ 2 1 t _ t 

It can be verified that in the coordinates (T, p, 6, (j)) the hne element of the near-horizon extremal black hole metric 
becomes 

ds^ = dT^ + dp^^ + C{r+f{d0^ + sin^ 0d4>'^), (20) 

which is a Bertotti- Robinson type of the spacetime. Making use of the near-horizon approximation we obtain the 
following value of the spinor function: 

-0-1 1 r, (21) 



where 7 implies 



Now, the spinor function diverges near the black hole event horizon as (r — r+)'''~^/^, where 7 depends on the metric 
coefficients taken at r+ . In the case under consideration, one has the influence of a global monopole parameter b on 
the divergence of the zero mode spinor function. In the case of the extremal dilaton black hole the situation is much 
more complicated. As can be seen, for the extremal dilaton black hole the event horizon of it is singular in Einstein 
frame and has vanishing area. Just C(r+) — )• in this limit. However, in string frame previously singular horizon has 
been pushed off to an infinite proper distance. The very similar situation was indicated in the case of the expulsion 
of the Higgs vortex from an extremal dilaton black hole [36| . In the case under consideration a full numerical study 
would be required. To sum up, the divergence of spinor function for zero modes occurs in the near horizon limit as 
was expected from the point of view of the no-hair theorem. 

Returning to relations ([TT |) - p^ . we observe that they can be rewritten in the form of a coupled second order 
differential equations system. Consequently, one obtains 

dU-dlj + H,{r)f = Q, (23) 
d^,g~dlg + H2{r)g = Q, (24) 

where we have denoted by Hi{r) = H'^{r) -t- dr,H{r) and H2{r) = H^{r) — dr,H{r) the effective potentials for the 
s-wave sector while the potential H(r) = Mr)(.r-a{r)c{r)) ^ 

One can remark that Hi{r) and H2(r) are supersymmetric partners, in the sense presented in Ref. [40j . derived from 
the same superpotential H{r). In Ref.[4l| it was proved that these relations between potentials provided that they 
were the sources of the same spectra of quasinormal modes. It means physically that Dirac fermions and antifermons 
have the same quasinormal modes in the considered general spherically symmetric black hole background. 



B. fc > modes for Dirac fermions 



It happened that the eigenspaces with eigenvalues k and m are four-dimensional ones for fc > 1 j^. Operators 
appearing in the Dirac equation can be expressed as matrices in the basis (/l m,j+>, gL\k,m,— >, fji\k,m, j_ > 
, gnlk, m, + >) in which the operators in question are diagonalized. On this account, it can be readily seen that 

(fL\k,m,j+ >\ 
gB.\k,m,+ > I 

while the transverse Dirac operator implies the relation 

fc + 1\ 

(26) 
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Also, we may note that the matrix form of the operator na x t is given by the following: 



_ 2i ^/ki 

n<7 X T 



2fc+ 1 







^k{k + l) 

-k y/k{k + l) 

k 



k + 1 



-^k{k + l)J 



It can be verified that in this situation the Dirac equations reduce to the form 

dtfR + dr*gR + l3{r)gR + a{r)fL^0, 

dtgR + dr'fn - /3(r)/fl - a{r)gL = 0, 

dth - dr-gL - P{r)gL ~ a{r)fR = 0, 

dtgL - dr'fL+ P{r)fL+a{r)gi^^^O, 

where by a{r) and /3(r) we have denoted the following quantities: 

Mr) _ a{r) A{r) 



C{t) 



r{2k + 1) 



(27) 



(28) 
(29) 
(30) 
(31) 



(32) 
(33) 



In Ea. (|55| Dk is equal to fc + 1 for the i?- function and equals to k for the L- functions. In order to simplify the above 
equations we put for all the functions an explicit time-dependence in the form exp(—iujt). We assume further, that 
the following relations between considered functions = ifn and gj^ = —ign are fulfilled. It provides the following: 



d 

dr* 



A{r) 



-l3{r) ia 
ia l3{r) 



0, 



where a(r) and /3{r) are given by 



A{r) 2k + 1 

r {2k + 1) - a{r) C{r) 
(2A: + 1) C{r) 

In order to simplify further the radial Eq. (p4l) let us make a change of variables defined as 



= T 



sm ■ 



(34) 

~t N a{r) ^k{k + l) 

"(0 = THT—, ' (35) 

/3('^) = ^^^^^^TT^Vt^TiP^ (^ + 1)- (36) 

(37) 

arctan ( — | . (38) 

(39) 

It happens that they may be simplified further. Namely, let us make another change of variables given by dr = 
Consequently, we arrive at 



where the angle 6 yields the following: 



/ ia 



Hence, the underlying relations reduce to the form as 



d 

dr^, 



A{r) 



/32 



-1 



1 d0 

2uj dr.. 



uj 



a 



W{r) 







df \f 

where by W{r) we have denoted the following expression: 

2ujA{r 



-cj 
UJ 



W{r) 



(10 



(40) 



(41) 
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The set of Eas. (p9)) can be decoupled providing the the system of second order differential equations. Namely, they 
yield 

- Wi{r) -g^L,' ~g, (42) 
~ W2{r) f = -uj' f, (43) 

where the effective potentials Wi (r) and W2 {r) imply 

Wi{r) = W^(r) + dfW{r), (44) 
W2{r) ^W^{r) - dfW{r). (45) 

Because of the fact that Wi{r) and W2{r) are supersymmetric to each other, g and / will have the same spectra both 
for quasinormal modes and scattering. 

C. Massive Dirac Fermion Modes 

So far we have considered massless fermion case. Now, we wish to generalize our considerations and present some 
relevant arguments concerning with a massive Dirac fermion. The equation under consideration implies 

i 7''(Vp - iXH^)ij - = 0, (46) 

where m is a mass of the Dirac fermion. By the same procedure that we followed in the preceding section, the Dirac 
equations for the massive case can be written as 

dtfR + dr.gR + f3{r)gR + a(r)/i + i m A{r)fL = 0, (47) 

dtgR + dr^fR - Pir)fR - a{r)gL + i m A{r)gL = 0, (48) 

dtfL - dr^gL-m9L-a{r)fR + imA{r)fR^0, (49) 

dtgL - dr^fL+P{r)fL+a{r)gR + im A{r)gR=0, (50) 

Further, assuming that /l — ifR and gL = —ign., we obtain 



dr*\gj I ai -mjKg) \mAir) 1 [g 



(51) 



where ai(r) is provided by 



r(2/c + 1) 

As in the previous section, let us introduce 9(r) defined as 



. a{r) A{r) ^k{k + l) , 

Qfi = i[uj — — . (52) 



{/3(r)) 



9{r) = arctan[ ^ ) . (53) 

Applying the transformation T and introducing new variable defined as dr — ( A{r) — (l/2m)^ ) dr^, we arrive at 
the following expression: 



± ff\ _ v//3(r)^+ai(r)2 f / 

j^fy.) L d9_ 



dr \gj f 1 de \ \-gJ \-f 



g 



(54) 



It can be also shown that these equations decouple. Namely, one gets 



dj f - Gi(r) f + m^ / = 0, (55) 
dl~g - G2ir) g + m^ g = 0, (56) 
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where Gi(r) and G2{r) yield 

Gi{r) = df.W + W^, G2{,r)^^dfW + W^, (57) 

while by W{r) we have denoted 



^- A{r)-^^ ■ ^^^^ 

The above potentials have the same features as potentials in preceding sections. Namely, they are supersymmetric to 
each other. 



III. THE DECAY OF DIRAC FERMION HAIR 



This section will be devoted to the problem of the decay of massive Dirac fermion hair in the background of 
black holes with non-trivial topology of the event horizon. We shall use the sign convention presented in Ref. [27). 
where it was pointed out that that the treatment of fermions in spherically symmetric backgrounds may be greatly 
simplified by recalling a few basic properties of the Dirac equation. Namely, for a line element of the form as 
g^^dLx'^dx" — gab{x)dx°'dx^ + gmn{v)dy'^ dy"- the Dirac operator.^ = 7^V^ can be decomposed as a direct sum 



X • -r y 



(59) 

Moreover, under a Weyl conformal rescaling of metric tensor g^i, = ^^g^n it could be shown that 

pij = i7-i("+i)|)v; , %p = (60) 

Let us consider a conformo-static metric of the form as 

ds^ = -A^df + ^^dx'dx' , (61) 
where A = A{x^) and $ = ^{x/), i ~ 1,2 . . .n ~ 1. Then, one obtains 



and finally we may note that 



ds' = A^ ( -dt^ + {^^dx'dx' ] , (62) 



where^ is the Dirac operator of the metric (^^^^dx^dx^ and '0 = yls^" ^'^ip. Using again the conformal property one 
arrives at the following: 

n^{^y-^''''\'dd, (64) 

with -ip — (^) ^ ^^■(/i. Let us suppose moreover, that ^I^ is a spinor eigenfunction on the (n — 2)-dimensional transverse 
manifold VL satisfying relation 

p,,^ = A*. (65) 
In the case of a (n — 2)-dimensional sphere, the eigenvalues of the spinor ^ were given in Ref. in the form 

A^=(/ + ^)', (66) 

where Z = 0, 1, . . . 

Furthermore, one can assume that p^j = mtp and set the following form of ip: 
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It can be checked by the direct calculations that the above form of fermion fields provides the following: 



(68) 



One should recall that the matrices 7", 7^ satisfy the Clifford algebra in two spacetime dimensions. If we assume that 
■0 oc e~"^* it can be shown that we are left with the second order equation of the form 

0+c.2x^A(.)2(--^)^. (69) 

Although the above derivations are valid for arbitrary number of spacetime dimensions, in what follows we shall 
restrict our attention to the four-dimensional case. 



A. The background of a black hole with a global monopole 

We first focus on the case of decaying massive fermion hair on the background of a black hole with a global monopole. 
In four-dimensional spacetime the line element describing a black hole which swallowed a global monopole is written 
as 

ds^ = - fl - SirGif - ^) + — ^ + r^dn\ (70) 

where by M we denote mass of black hole and rj is the symmetry breaking scale when the monopole is produced. If 
we introduce the coordinate transformation in the form 

< ^ (1 - 87rG7?2)-3/2t^ r ^ [l-^TTGr^^y^'^r, (71) 

as well as the new parameters which yield 

M = {l-^TTGif Y^I^M, b"^ = I - SnGr)^ , (72) 

then, the line element of a black hole with global monopole reduces to the form given by Eq.(IT]). In what follows we 
put G = 1 for simplicity. 

The spectral decomposition method will be our main tool in the analysis of the time evolution of a massive Dirac 
spinor field in the background of a black hole with global monopole. It was revealed in Refs.ji^ that the asymptotic 
tail was connected with the existence of a branch cut situated along the interval —m < lu < m. An oscillatory inverse 
power-law behaviour of a massive Dirac field arises from the integral of the Green function G(j/,?/';w) around the 
branch cut. The time evolution of the massive Dirac field takes the form 



Xiy,t)^ j dy' G{y,y';t)xt{y',0) + Gt{y,y';t)x{y',0) 

for t > 0, where the Green's function G{y,y';t) implies the following relation: 

92 92 



(73) 



, V 



G{y,y';t)=S{t)S{y-y'), (74) 



where V is an effective potential. 

In order to find the Green function in the case under consideration we shall use the Fourier transform (43j G{y, y'',u}) = 
J^dt G{y,y';t)e^'^* and reduce Eq.dZll) to an ordinary differential one. The Fourier's transform is well defined for 
Im w > 0, while the inverse transform becomes 

-1 />00+Z€ 

G(2/, y'; t) = — doj G{y, y'; w)e-*-*, (75) 

'^'^ J-oo+ie 

for some positive number e. Thus, in this picture the above Fourier's component of the Green function G{y, y'] lu) can 
be rewritten in terms of two linearly independent solutions of the homogeneous equation of the form 

-^+u;^-v)x.^O, « = 1,2, (76) 
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where we have denoted hy V = A{r) yni — 
Let us consider the boundary conditions of the problem in question. For Xi they are described by purely ingoing waves 



crossing the outer horizon of the static charged black hole xi 



^'^y as y — ?> — oo. As far as X2 is concerned, it 
should be damped exponentially at i+, i.e., X2 — e"^™ ~" ^ at y — > oo. 

In our considerations we shall assume that the observer and the initial data are situated far away from the black 
hole with global monopole. Eq. (j76|) may be rewritten by using new variables expressed as 



1 - 



2M 



1/2 ' 



(77) 



where i = 1,2. Next, we expand Ea. ([7S)) in a power series of 2M/r neglecting terms of order ©((w/r)^) and higher. 
Having this in mind we are left with the equation 



2 2 

— m + 



4A/6a;2 + 2Am 
br 



2 1 



^2y,2 



(78) 



The main result of the above procedure is the conclusion that Ea. ([75)) may be solved in terms of Whittaker's functions. 
Two basic solutions are needed to construct the Green function, with the condition that \ u \> m. The Whittaker's 
functions provided the solution of the above equation are xi = Ms^fi{2ujr) and X2 = Ws.jx{2LL:r), while their parameters 
are given by 



1/4- 



A2 



S = 



iMbuj^ + 2 Am 



52 ' " 2w6 
It is sufficient to conclude that the spectral Green function takes the form as 



~ 2 2 2 

Ul =171 — U! . 



Gc(r,r';i) 



1 

2^ 
1 

2^ 



dw 



Xi(r,J)e''*) X2(?-',<ie") Xi{r,uj) X2{r',uj) 



dwf{Cj) e-™*. 



W{Cj) 



(79) 



(80) 



where W{Cj) is the Wronskian. 

We first focus our attention on the intermediate asymptotic decay of the massive Dirac hair, i.e., in the range of 
parameters Af ^ r <C t <C M/{mM)'^. The intermediate asymptotic contribution to the spectral Green function 
integral gives the frequency equal to a) = 0{^m/t), which in turns implies that (5^1. Having in mind that 5 
results from the 1/r term in the massive Dirac field equation of motion, it depicts the effect of backscattering off the 
spacetime curvature and in the case under consideration the backscattering is negligible. Taking into account all the 
above we obtain the result 



22A-ir(-2/i) r(i + /2) 
Ar(2A) r(i - A) 



J2A+l)7ri 



(rr ) 



where we have used the fact that dir ^ 1 and the form of /(w) can be approximated by means of the fact that 
M{a, 6, z) = 1 as z tends to zero. The resulting Green function reduces to the form as 



23'i-l r(-2/i) r(i + ^)r(^ + i 



1 



„(2A+l)7ri 



(rr') - 



flT{2il) r(i - m) 

In the limit when t 3> 1/m one can show that the spectral Green function yields 

23A-ir(-2/i) r(i + A)r(/i + i) 



-^i+u("ii)- 



(82) 



Gcir,r';t) 



fi^ ilT{2fi) r(i - m) 



1 + e(2'i+i)'^' ) (rr') m'^ r'^-'' cos(mi - -{jl + 1)). (83) 



Eq. ()83p depicts the oscillatory inverse power-law behaviour. In our case the intermediate times of the power-law tail 
depend only on fl which in turn is a function of the multipole number of the wave modes and monopole parameter b. 
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On the other hand, for the late-time behaviour quite different pattern of the decay should be expected because of 
the fact that backscattering off the curvature of the spacetime play an important role. For this case k ^ 1 and /(w) 
may be rewritten using the fact that Ms^p,{2ojr) w r(l + 2/i)(2wr)i (5^'^ J^{y/8Suir). It yields 



r(i + 2/i) r(i - 2/2) 

2/i 

2Ar(2/i) T{^-~^^5) 



(84) 



+ e(2/''+i)/2^(V8fc^) /2^(\/8fc7) 



The first part of the above Eq. (|M| the late time tail is proportional to t ^ and it occurs that we shall concentrate on 
the second term of the right-hand side of Ea. (|84)) . For the case when k ^ 1 it can be brought to the standard form 



written as 



where we have defined 

while N provides the relation as follows: 
(r(l + 2A))2 T{-2fi) 



1 + (-l)2Ae-2^«'5 



TV : 



2/i T{2jl) 



(rr')^ 



J2a(V8^) J2~^{^m^) + hf^iVmJ^) /2^(V8fcr') 



(85) 



(86) 



(87) 



At very late time both terms e™* and e^'"^ are rapidly oscillating. It means that the scalar waves are mixed states 
consisting of the states with multipole phases backscattered by spacetime curvature. Most of them cancel with each 
others which have the inverse phase. In such a case, one can find the value of Gc(2) by means of the saddle point 
method. It can be found that the saddle point is given by 



ao 



TT (4M6a;2 -|- 2\m) 



2V2bT 



In comparison to the late-time behaviour of the second term in Eq . ((84|) . the first term can be neglected. The dominant 
role plays the behaviour of the second term, i.e., the late-time behaviour is proportional to — |. Consequently, this 
fact implies the resultant form of the spectral Green function for the late-time behavior of massive Dirac field in black 
hole with global monopole spacetime, namely it can be written as 



Gc(.,/;t) = ^(.)t 



AMbm^ + 2Xni 



(mt) 6 sin(TOi) x(r', m) x(r', m). 



(89) 



One can see that the late-time behaviour is independent of the global monopole parameters as well as the mass of the 
Dirac field. 

B. The background of a black hole with a cosmic string 

In this subsection we analyze the decay pattern of massive Dirac hair on the spherically symmetric solution of 
dilaton gravity being the low-energy limit of the string theory. In four spacetime dimensions, the action for the 
dilaton gravity with arbitrary coupling constant a is given by 



S — I ct^x^J—g 



(90) 



where (/> is the dilaton field, = 2V[^A^] is the strength of the U{1) gauge field. 

This theory constitutes the low-energy approximation of the heterotic string theory. The static spherically symmetric 
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black hole solution with a cosmic string passing through is a thin string approximation of the Nielsen- Olesen vortex 
piercing the black hole. The metric for such a system was found in Ref.|30|. It was assumed that the system underlies 
the complete separation of the degrees of freedom between the U{1) gauge field responsible for the charge of the black 
hole in question and the gauge vortex field. 

The line element we shall consider will be given by Eq.([2]). The same procedure as in the preceding section and 
Ref.(27l| enables us to find the intermediate and late-time asymptotic behaviour of massive Dirac hair in the background 
of dilaton black hole with a cosmic string. For the completeness and the reader's convenience we describe the crucial 
points of the underlying procedure. Namely, it is convenient to change variables in Eg.lfTS]) in the way as follows: 

X^ = -T^, (91) 



1-^ 1- 



2(1 + q2) 



where i — 1,2. Neglecting the terms of order 0{{uj/r) ) and higher leads us to the relation for ^ 



w — m H (92) 

r 



X — 2Amr_(ai + 02) + m air^r 



which in turns provides the solutions given by Whittaker's functions. Namely, xi = Ms^jx{2ujr) and X2 = Ws^fi{2ujr) 
with the following parameters: 



= yjl /4 + A2 - 2Xmr- + m?air+r^ , (93) 



5 = 



a;2(r+ + air_) + Am(l + r+) - ^(r+ + air_) 



- 2 2 2 

cj — m — u] . 



The intermediate asymptotic behaviour of the massive Dirac hair on the dilaton black hole pierced by a cosmic string 
will be given by Eq. (j83p . but in this case the parameters of Whittaker's functions fi and 5 are of the form provided 
by the relation (|93|) . 

As far as the late-time asymptotic behaviour is concerned it yields 



G,(r,/;i) = ^m2/3 (^)l 



2m^(r_|_+air_)+2Am(l+r+)— r7i^(r++air_) 



{mt) 6 sin(mt) %(r, m) x(^'j 'Ti). (94) 



It envisages the fact that the late-time asymptotic decay pattern of massive Dirac hair in the background of spherically 
symmetric dilaton black hole is proportional to —5/6. 

Having in mind the properties of the Dirac operator presented in the preceding sections, we can see from the 
relations (|65p - (l67p that the main modification will stem from the quite different topology of the transverse manifold. 
In our case it will be S'^-sphere with a deficit angle. The reason of it is the presence of a cosmic string passing through 
it. The other part of the Dirac function will be not affected by the string. Now we proceed to the analysis of the 
eigenvalues of the Dirac operator on the underlying manifold. 
The metric on the sphere with cosmic string has the form as; 

dn^ = B^sm^0(f(f) + (f9, (95) 

where B = 1 — 4/z, while /x is the mass per unit length of the cosmic string [s^ . 

As can be seen by the direct calculations the Dirac operator may be written as follows: 

^^,V = 7'a, + -^9^^ + ^^^[7^7V, (96) 
B sm 4 



where gamma matrices satisfy the Clifford algebra rules in two-spacetime dimensions. 



The square oipl^2 takes the form 



2 (9e(sin6'ae) , 1^2 cos6' . 11 
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Let us assume that the eigenfunctions of the Dirac operator (|96l) are two-component spinors that fulfilled the following: 



Me, 4 



E 



V2 \P>^m{e) 



(98) 



where m are half-integers. On evaluating the Dirac equation on a stringy sphere we find that it provides the system 
of differential equations expressed as 



de + 
d0 + 



cot ( 
cot t 



a.Xm{0) 



m 
sm6 

m 
sin 6 



where rh = m/B. Next, we change the variables x = cos 
I ^ 2 ^ '"■ ~ '"- '-'3 ^ 

dx 



It allows us to rewrite the underlying relations as 



1 



f3xm{0) 



(99) 
(100) 

(101) 



where is the Pauli matrix. The above equations are singular at the poles of the sphere .t = ±1, so it is convenient 
to redefine the unknowns 



axm{d) 

Pxm{0) 



(l_a;)5lA-5l (i + a;)5l™+5l ixm{x) 
(l-a;)5l™+5l (l + a;)5l*-^l 



■nxmix) , 

It can be verified by evaluating these expressions, that one can achieve to the relations 



l-x' 



dx^ 



CT3 



2 I m I 4-2 ) x] — - m{m + 1) + { \^ - - 
' ' dx \ A 



VX,n[X) 



(102) 



(103) 



Having in mind the general formula for Jacobi polynomials (see e.g., |44| ) one gets the square of the eigenvalues for 
the Dirac operator on the S'^-sphere with a cosmic string passing through it, which yields 



X'=[l+\ fn I +i 



(104) 



where I — 0, 1, 2 . . . Because of the fact that fh ~ m/B the crucial role is played by the factor B connected with a 
mass per unit length of a string. 

To conclude this section we remark that recently the non-linear origin of the power law tail in the long-time 
evolution of a spherically symmetric self-gravitating massless scalar field was discussed ([i^ and references therein). 
The perturbation method was used to obtain the expression for the tail and then numerical integration was performed 
to check the results. The non-linear evolution of a black ringdown in the framework of the higher-order metric 
perturbation theory was conducted in Ref. [46| . It was argued that these non-linear components should be particularly 
significant for binary black hole coalescences. 

Although, it is undoubtedly that the linear theory is clearly useful [i^ H^l, the recent numerical simulations of 
the aforementioned problems yield that non-linearity is also worth elaborating. We hope to return to the problem of 
non-linear behaviour of massive Dirac hair elsewhere. 



IV. CONCLUSIONS 



To summarize, we have analyzed the behaviour of Dirac fermions in the background of non-trivial topologies. 
Assuming the complete separation of the degrees of freedom of the fields in question, we have first considered the 
massless zero modes in the near horizon limit for both nonextremal and extremal black holes. We have found that the 
spinor function is divergent near horizon as (r — r+)^^/^ for the non-extremal case. In the vicinity of extremal black 
hole horizon one gets the dependence (r — r+)'''^^'*, where 7 was built from adequate components of line element 
taken at r+. It was also established that the global monopole parameter influenced the divergent behaviour. Namely, 
the bigger it is the more divergent the spinor wave function is. On the other hand, in dilaton gravity black hole 
spacetime one has the coupling constant influence on the divergence of the spinor wave function. 



14 



Then, we take into account Dirac fcrmion modes for k > and show that the Dirac equations can be decoupled to 
the system of second order differential equations. We treat also the massive case of Dirac fermion fields. As in the 
massless case one has also the situation that underlying equations decouple to the system of second order differential 
equations with the so-called supersymmetric potentials. 

Next, we proceed to the intermediate and late-time behaviour of massive Dirac hair in the backgrounds of black 
holes with topological defects. Namely, we considered the black hole with global monopolc and cosmic string passing 
through it. In the case of black hole which swallowed a global monopole one gets the modification of the intermediate 
late-time behaviour which depends on mass of the field in question as well as the global monopole mass. The 
intermediate late-time decay of the hair is quicker comparing to the decay rate of massive Dirac hair on a black hole 
without global monopole. The intermediate oscillatory power-law depends also on the multipole number of the wave 
mode. But it is not the final pattern of the decay of the adequate massive hair. At very late times the resonance 
backscattcring off the spacetime curvature plays the dominant role. This decay pattern is independent of the presence 
of a monopole and of the mass of the Dirac field, and it is proportional to t~^^^. 

The analysis of the decay of massive Dirac hair on the background of a dilaton black hole with a cosmic string 
passing through it, reveals that the intermediate as well as the late-time behaviour is independent on the presence of 
this kind of topological defect. The main modification appears in the eigenvalues of the Dirac operator on a S^-sphere 
pierced by the cosmic string. The cosmic string parameter connected with its mass per unit length plays the crucial 
role in this case. On the other hand, the intermediate late-time behaviour depends on the multipole number of the 
wave mode as well as the mass of the Dirac field. In turn, the late-time behaviour of the massive Dirac field is 
independent of the above factors and it is proportional to t"^/^. 

Having in mind the previous works treated the scalar, fermion and vector black hole hair decays one can conclude 
that at asymptotically late-times the decay of the hair in question is universal and does not depend on spin of the 
field, wave number of the mode and the topology of the black hole event horizon. 
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